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Abstract 



We studied in this work a modified version of tlie Ziff, Gulari and Barshad 
(ZGB) model in wliich the adsorption probabihties of CO and O2 are site 
dependent. We employed Monte Carlo simulations to determine the phase 
diagram of the model as a function of the partial pressure of CO in the gaseous 
phase, and of the relative adsorption probabilities. For the case where there 
is no site with affinity only to O2 adsorption, we found that the non physical 
0-poisoned state disappears. This feature can not be seen in the usual ZGB 
model. 

PACS: 05.70.Ln, 82.65. Jv 
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I. INTRODUCTION 



The search for models to explain the dynamic behavior of nonequilibrium physical sys- 
tems which present heterogeneous catalytic reactions are the aim of several works in statis- 
tical physics [1-9]. The well-known Ziff, Gulari and Barshad model (ZGB model) [10], used 
to explain the kinetics of reactions hke CO + 1/202^^02 and 2CO + 2N0 N2 + 2CO2 
[5,ll]on a catalytic surface, shows some features that are not in agreement with the exper- 
iments. For example, this model presents, for the catalytic oxidation of CO, a continuous 
phase transition between an active and an 0-poisoned states, which is not found experi- 
mentally. Several experiments tell us that the production of CO2 molecules starts as soon 
as the CO molecules begin to reach the surface [12]. Some attempts were made with the 
aim to clarify this point [2,5,6,9,12,13]. For instance, Hua and Ma [9] proposed a model 
where the sites on the catalytic surface are separated in two types: active and not active 
sites. In their model, the surface reactions follow the Langmuir-Hinshelwood mechanism, 
although the adsorption of molecules depends on the kind of the chosen site, i.e., if the 
site to do the adsorption is active or not active. They prepared the lattice surface with a 
fraction / of distributed randomly active sites. After, they considered that the conditions 
for the adsorption of molecules on the surface are: (a) an O2 molecule is adsorbed if at least 
one of the two vacant adjacent sites, randomly chosen, is an active site; (b) a, CO molecule 
is adsorbed with probability 1 if the empty site, randomly chosen too, is an active site; if 
the chosen site is not active, the adsorption of CO occurs with probability p. They found 
the phase diagram of their model and showed that the phase transition from the reactive 
state to the CO-poisoned one is a continuous transition for small values of /, and it is a 
discontinuous phase transition when the value of / is large. When they increase the value 
of p, the reactive region is narrowed. Besides, they showed that the O-poisoned phase does 
not appear for different values of /. 

In this work, we propose a variation of the ZGB model. We aimed to study the influence 
in the phase diagram of the hipothcsis intoduced by Ziff et al. that all sites at the catalytic 
surface are equally accessible to adsorption for both incident molecules {CO or O2). In our 
model, the catalytic surface does not present a site distribution with equal probability for 
adsorption for both molecules: we proposed that the adsorption of a CO molecule only occurs 
if the chosen empty site has preference for the adsorption of this molecule, and analogous 
procedure is followed for the adsorption of O2 molecule. For this, the lattice presents a 
fraction Aco of sites that prefers adsorb CO molecules, a fraction Aq of sites that prefers 
adsorb O2 molecules, and a fraction Af, of sites that accepts both molecules. This model can 
be understood as a catalytic surface with deffects, where the affinity to adsorption is site 
dependent. 

The model was studied by Monte Carlo simulations, and we determined the production 
rate of CO2 as a function of Aco, Aq and of the relative adsorption rate of CO molecules, 
denoted by yco- In this problem, the significant parameter to tell us if the system is in 
a reactive state is the production rate of CO2, Pco-z, because a nonzero number of empty 
sites does not guarantee that the system is active - it is necessary to know what kind of 
molecule the vacant sites at the lattice are available to adsorption. A possible situation to 
be found is that where we have only one empty site in the lattice, and this site accepts only 
the adsorption of oxygen: once the O2 molecule needs two vacant sites to be adsorbed, this 
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site always will be empty. With the purpose to rich our model, we separate the sites in the 
lattice in three kinds: a fraction Aqo of sites that accepts only the adsorption of CO, a 
fraction Aq of sites that accepts only the adsorption of oxygen, and a fraction of sites 
that accepts the adsorption of both molecules. The ZGB model would be a particular case 
of our model, in which the fraction Ai, would be 1. We show that the present model predicts 
a phase diagram where an absorbing phase by oxygen atoms does not appear when Aq = 0.0 
and Aco varies from 0.0 to 0.8. In the next Section we present the model, and in Sec. Ill 
we present the results obtained for the production rate of CO2, the discussion of the results 
and conclusions. 

II. MODEL AND MONTE CARLO SIMULATIONS 

The simplest model to describe the surface catalytic reaction is that one presented by 
Ziff, Gulari and Barshad [10]. In this model, the catalytic surface is represented by a square 
lattice, and all sites are accessible to molecules of any kind. In this surface, molecules of 
two species are adsorbed: CO and O2 molecules. When a CO molecule is adsorbed, it 
occupies only one site in the lattice, while the O2 molecule occupies two sites because it 
dissociates in two oxygen atoms when it is adsorbed. If the adsorbed CO has an O atom 
as its nearest-neighbor, CO and O react forming a CO2 molecule. This new molecule leaves 
immediately the surface, and two sites are set free in the lattice. The whole process follows 
the Langmuir-Hinshelwood mechanism and the following three steps are considered: 

(1) CO(g)+V^CO(a), 

(2) 02(g) + 2V^ 20(a). 

(3) CO(a) + 0(a) ^ C02(g) + 2V, 

where the labels a and g denote adsorbed and gaseous particles, respectively, and V is a 
vacant site. The steps (1) and (2) describe the adsorption of the molecules CO and O2 , 
respectively, and the third step represents the proper reaction between the adsorbed species 
to form the CO2 molecule. In this model, the vacant sites are free to adsorb both molecules 
equally. When we take into account that the sites in the lattice are divided into three 
categories, those steps may be changed: 

(1) CO(g)+Vco,b^CO(a), 

(2) 02(g) + 2Vo,b^ 20(a), 

(3) CO(a) + 0(a) ^ C02(g) + Vco,bVo,b, 

where the labels CO and O on V indicates if the vacant site K is a preferential site for 
a CO or O2 adsorption, respectively. The label b means no preferential adsorption site. 
To describe the whole process we need to consider also the relative adsorption rate of CO 
molecules, denoted by yco- 

The Monte Carlo simulation was performed on a square lattice of hnear dimension L. For 
a given value of L, the lattice was prepared taking into account the fractions Aco, and Aj, 
of sites that accept the adsorption of the molecules. We have used in all of our simulations 
periodic boundary conditions. For each value of yco, we choose what kind of molecule will 
be tried to adsorb. If the chosen molecule is CO, we choose randomly a site a where CO 
can be adsorbed: if this site is empty and if it accepts to adsorb this molecule, i.e., this is 
a vacant site of the type Vco or V5, the adsorption is performed; if a is occupied, the trial 
ends. If the selected molecule is O2, we select randomly two nearest neighbor sites (7j and 
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af if both sites are empty and if they are of types Vq or the adsorption happens; if one 
of them is occupied, or one of them is of Vco type, the trial finishes. Since the adsorption 
of CO occurs, its four nearest neighbors are analyzed: if one of them is an O atom, a CO2 
molecule is formed, and it leaves immediately the surface. When this happens, two vacant 
sites are free for a new adsorption. The same occurs when an O2 molecule is adsorbed: in 
this case can be formed one or two CO2 molecules, with two or four sites being free in the 
lattice, respectively. One Monte Carlo step (MCS) equals L x L adsorption trials. After 
one MCS, we compute the production number of CO2. In order to estimate the quantity of 
interest, we discarded the first 1000 Monte Carlo steps to reach the stationary state, and we 
have used the next 6000 Monte Carlo steps to calculate the averages. 

In our model, the significant parameter to determine if system is in a reactive or in a 
non reactive state is the CO2 yielding, Pco2 (iii the ZBG model, we can take the fraction 
of the vacant sites in the lattice). The parameter Pco2 used because, even having a 
nonzero fraction of vacant sites in the lattice, we can not guarantee the state is active: it is 
necessary to know for what kind of molecules the vacant sites are available to adsorption. 
For example, if in the lattice there is only one vacant site and it does not accept CO 
molecules for adsorption, the site will never be occupied. Figure 1 shows some snapshots of 
a square lattice with L = 50, for yco = 0.45, Aco = = 0.3 and Ah = 0.4, for the inicial 
configuration (Fig. la), after 100 MCS (Fig. lb) and after 1000 MCS (Fig. l.c): we can see 
that the system always presents a number of vacant sites {Vco, Vo a-nd Vb) different from 
zero. 

III. RESULTS AND CONCLUSIONS 

In the ZGB model, the phase diagram shows two absorbing phases (poisoned or inert) 
and one active (or reactive) phase. The transition from the 0-poisoned phase to the reactive 
one occurs for yco — 0.389, and this is a continuous phase transition. If the value of yco 
increases, the system remains in a reactive phase up to yco — 0.527, when it becomes 
CO-poisoned. This second transition is a first-order one, where there is an abrupt change 
in the number of vacant sites for yco = 0.527. However, when the lattice sites present 
different affinities for the adsorption of molecules, the phase diagram changes. If Aco or 
Ao are different from zero, the inert phases {CO or O-poisoned phases) disappear. In the 
sequence, we present some representative diagrams of the simulations, all the results for a 
square lattice with L = 50. 

Figure 2 shows the results for the production rate of CO2, Pc02i f^o^ ^co = 0.0 and Aq 
from 0.0 to 0.8 (we remember that Aco + Ao + A^ = 1). We can see that Pco2 decreases 
drastically when the value of Ao increases, becoming almost zero when Ao is close to 0.8. 
The figure shows that the maximum production of CO2 occurs at different values of yco, 
and its value decreases when yco increases; also, when Ac is increased, this maximum 
moves in the direction of increasing values of yco- In the case present, the CO-poisoned 
phase does not occur. Figure 3 exhibits the production rate of CO2 as a function of yco, 
for Ao — 0.0 and Aco assuming values from 0.0 (ZGB model) to 0.8. This is the most 
interesting case because it is related to disappearance of the non-physical O-poisoned phase, 
which is predicted by the ZGB model. Also, the character of the transition from the reactive 
phase to the CO-poisoned one changes: the first-order phase transition becomes a continuous 
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transition when the value of Aco is larger than 0.2. Contrary that observed in Figure 2, 
the maximum of Pco2 assumes increasing values when yco also increases. An intermediate 
case is that shown in Figure 4: it displays Pco2 ^ ^ function of yco, for Aq = 0.5 and Aqo 
assuming values in the range 0.1 to 0.5. We can see, in this case, that the system does not 
present any inert regions, and the maximum production of CO2 is almost constant, only 
changing the value of yco, for which the maximum occurs. 

In summary, we have shown that with a slight modification of the ZGB model, which 
includes a non equal probability of adsorption of the species, it is possible to account for the 
experimental result concerning the absence of 0-poisoned states in the oxidation of CO. We 
have found the production rate of CO2 as a function of the relative adsorption rate of CO 
molecules for several conditions of the catalytic surface: the probability of adsorption is site 
dependent and this fact leads to changes in the phase diagram. The most interesting case 
occurs when there is no lattice site with affinity only to O2: the non-physical O2 poisoned 
state, predicted by the ZGB model, disappears in the present model. We stress that the 
model can be understood taking into account deffects that modify the chemical properties 
of the lattice sites. 
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FIGURES 



FIG. 1. Snapshots in a square lattice oi L = 50, and for Aco = Aq = 0.3. (a) Initial 

configuration. Legend: Up triangle corresponds to the sites that are accessible to the adsorption 
of O, Dark down triangle correponds to the sites that accept only CO molecules, and blank spaces 
accept both molecules, (b) Lattice after 100 MCS. Legend: Asterisks corresponds to a CO adsorbed 
molecule, Letter O to a O adsorbed atom and blank space corresponds to a vacant site, (c) Lattice 
after 1000 MCS. Legend is the same to (b). 

FIG. 2. Production rate of CO2 as a function of yco, for ^CO = 0-0 and Aq varyng from 0.0 
to 0.8. The lines serve as a guide to the eyes. 

FIG. 3. Production rate of CO2 as a function of ycOj for -^O = 0-0 and Aqo varyng from 0.0 
to 0.8. The lines serve as a guide to the eyes. 

FIG. 4. Production rate of CO2 as a function of yco, for = 0.5 and Aco varyng from 0.0 
to 0.5. The lines serve as a guide to the eyes. 
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FIG. 1 (a)-FULCO + GRANDI 
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FIG. 1(c)-FULC0 + GRANDI 
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FIGURE 2 



FULCO + GRANDI 
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FIGURE 3 



FULCO + GRANDI 




FIGURE 4 



FULCO + GRANDI 



